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Abstract: We investigate optical torques over absorbent negative refractive 
index spherical scatterers under the influence of linear and circularly pola- 
rized TEM 00 focused Gaussian beams, in the framework of the generalized 
Lorenz-Mie theory with the integral localized approximation. The funda- 
mental differences between optical torques due to spin angular momentum 
transfer in positive and negative refractive index optical trapping are out- 
lined, revealing the effect of the Mie scattering coefficients in one of the 
most fundamental properties in optical trapping systems. 
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1. Introduction 

Since 2000, when the first papers began to appear treating the subject of constructing some 
artificial medium with simultaneous negative permittivity and permeability [1,2], there has 
been an increasing interest on the new properties and revolutionary potential applications of 
what has been called negative refractive index (NRI) or double-negative (DNG) metamate- 
rials, left-handed (LH) materials or Veselago's medium (VM), which are artificial structural- 
arranged materials capable of delivering a homogeneous medium with an effective negative 
refractive index [3-5]. Some of these applications overcome current positive refractive index 
(PRI) limitations and, together with plasmonic structures, they are promising near-future 
technologies and devices for both microwaves and optics, such as in lenses, transmission 
lines, antennas, optical cloaking, cancer treatment and so on [6-14]. 

Recently, we have proposed the use of NRI metamaterials in optical trapping systems not 
only as optical devices for mechanical and lasing purposes, but as real trappable micro- or 
nano-particles. We have called this a "double-negative optical trapping" [15,16] or, alterna- 
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tively, "negative refractive index optical trapping" [17]. Radial and axial radiation pressure 
forces were then calculated for both Gaussian and Bessel beams using first a ray optics ap- 
proach, and further the generalized Lorenz-Mie theory (GLMT) with the integral localized 
approximation (ILA), thus allowing an all-optical regime analysis [16-19]. New and interest- 
ing trapping characteristics which could never be observed for any PRI particle were revealed 
for homogeneous and lossless simple NRI spheres. 

This paper extends our previous works and shows how the polarization state (linear or cir- 
cular) of a TEMoo focused Gaussian beam affects the optical torque exerted on lossy (absor- 
bent) NRI spherical particles in comparison to conventional PRI particles of the same geome- 
try. This kind of optical torque arises from the spin angular momentum (SAM) of the incident 
photons, in contrast to the torque produced by orbital angular momentum (OAM) transfer due 
to azimuthally asymmetric intensity profiles of incident beams such as Laguerre-Gaussian or 
high-order Bessel beams [20-22]. Due to the new resonances observed on the Mie scattering 
coefficients, every optical property (scattered and internal fields, forces, torques, scattering, 
extinction and absorbing cross sections and so on) is eventually affected in its amplitude or 
phase when the real part of the refractive index of the particle changes sign, the condition n rd 
= -1 playing no significant rule as does n rel = 1, where n rd is the relative complex refractive 
index between the particle and the surrounding medium. 

2. Optical Torques over NRI Spherical Particles 

An arbitrary electromagnetic wave can carry both spin and orbital angular momentum, the 
first being associated with its state of polarization, and the second with its azimuth light pat- 
tern dependence. The mechanism by which both are transferred to a given object is well es- 
tablished. It is well known, for example, that off-axis particles may rotate under the influence 
of a plane-polarized (linearly polarized) TEM 00 focused Gaussian beam due to an asymmetric 
linear momentum transfer [20,23]. Also, a circularly polarized laser beam carries an SAM of 
<7«/photon, where a = 0, + 1 or -1 for linear, right- and left-hand polarizations, respectively, 
and this AM causes a particle to rotate about its own axis, the sense of rotation being deter- 
mined by the state of polarization of the beam [20]. Finally, OAM of /w/photon is also carried 
by the fields, / being the topological charge or the mode index. A Laguerre-Gaussian beam, 
for example, has an AM of (cr + Z)«/photon [21,22]. 

Regardless of the kind of AM, Polaert et al. have shown that, in the framework of the 
GLMT, which is an extension of the Mie theory for plane waves [23,24], the Cartesian com- 
ponents of the optical torques exerted on a given object, whose centre coincides with the ori- 
gin of the rectangular coordinate system, are given by [23]: 
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where c is the speed of light in vacuum, M the refractive index of the medium surrounding the 
particle, k the wavenumber of the incident light, a n and b„ the Mie scattering coefficients 
(MSCs) for plane waves [25] and g™ m and g™ TE the TM and TE beam-shape coefficients 

(BSCs) for laser beams having an arbitrary spatial intensity distribution according to the usual 
GLMT nomenclature [26-29]. Reliable results demand an appropriate truncation of the 
double summation in Eqs. (l)-(3), thus demanding careful choices for p max and n max . One 
should note that all information regarding the impinging light is essentially encoded into the 
BSCs, whereas the geometric and electromagnetic properties of the scatterer are completely 
described by means of the MSCs. The expressions (l)-(3) are in accordance with previous 
theoretical works such as Barton et al [30] and are a slightly modified version of Eqs. (3)— (5) 
found in Ref. [23], where M stands for the refractive index of the particle probably due to a 
misprint. The direct use of Eqs. (3)-(5) from Ref. [23]. thus leads to incorrect torque predic- 
tions [31]. Expressions to compute the BSCs for plane waves and Gaussian beams with arbi- 
trary polarization using the ILA are found elsewhere and are not reproduced here [32]. A 
normalization factor of 2P/(nco 2 ), where P is the power and m the angular frequency of the 
incident wave, is omitted, in accordance with previous works [23]. 

Consider now that the scatterer is a homogeneous sphere with refractive index N = N re - 
iNi m , where N re can be either positive or negative depending on its PRI or NRI nature, respec- 
tively, and N im accounts for absorption. A time -dependence expO'co?) is implicitly assumed. 
The absorption causes the polarization of the beam to be changed after passing through the 
sphere, and SAM to be transferred to it [20]. 

In the next section, we show the behavior of the optical torque components for a plane- 
polarized (x -polarized) focused Gaussian beam with the particle being transversally displaced 
from the optical axis. Then, we observe the longitudinal z-component of the torque produced 
by circularly polarized beams due to a * 0, a situation that allows us to trap and rotate a NRI 
particle in a similar manner as currently performed with PRI microparticles. 

3. Numerical Results and Discussion 

3.1. Linear Polarization 

When a linearly polarized TEM 0 o laser beam hits a PRI homogeneous spherical particle, no 
torque is observed if the particle is located at the trap focus or, equivalently, at the point of 
stable equilibrium, regardless of its refractive index being real or complex. However, if the 
particle is transversally shifted along the trapping plane (perpendicular to the optical axis of 
the beam and containing the trap focus), it is well-known that an optical torque can be de- 
tected, due to an asymmetrical illumination, whenever this particle has a nonzero imaginary 
refractive index different from that of the external medium [20,23,30]. This N im 0 condition 
is also valid in the NRI case for achieving nonzero optical torques, as we shall see. 

To observe this angular momentum transfer, suppose an x-polarized Gaussian beam with 
wavelength A = 384 nm and beam waist radius a> 0 — 3.7 urn incident on an absorbent spherical 
particle with N im = 10~ 7 and radius a = 7.5 |im. The surrounding medium has M = 1 or, in oth- 
er words, N may be interpreted as the relative refractive index between the particle and the 
host medium. We assume that the permeability |i = 1 for a positive and |i = -1 for a negative 
refractive index particle. A first order Davis approximation is adopted [32-36]. Figures 1(a) 
and 1(b) show the radial torque component T x for several values of N r (> or < 0) and fixed 
longitudinal coordinate z 0 = 0 as the particle is shifted along y. These values for N r were cho- 
sen in order to facilitate the comparison with previous works with positive refractive index 
particles, especially when N r = + 1.19. In this case, T y is zero because of the relative position 
of the particle, whereas T, is always zero because of the polarization chosen [20,23,30]. Fig- 
ure 1(a) for N r = 1.19 with zo = 0 can be compared with Fig. 1 from Ref. [23]. Analogous re- 
sults were observed for T y when the particle is shifted along x. Figures 1(c) and 1(d) are simi- 
lar to the previous plots, except that now zo = 150 urn with different values ofN r . 
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y 0 (Mm) y 0 (urn) 

Fig. 1. T x profile for an A-polarized TEMoo laser beam (A = 384 nm, co a = 3.7 um) impinging on 
NRI or PRI particles with a = 7.5 um and \N im \ = 1CT 7 and different Nr. (a) and (b) (xoA>) = 
(0,0); (c) and (d) (x 0 ,Zo) = (0,150) um. In (c) and (d), the curves for \N,\ = 2.00 and 1.50 are al- 
most superposed and indistinguishable. The summation in (1) for these parameters leads to ap- 
proximate the same slopes for N, > 0 and N, < 0. The refractive index of the host medium has 
been normalized to M - 1 .00. 

By looking at Fig. 1, it can be inferred that, for the parameters chosen, the magnitude of T, 
for a NRI particle resemble that of the equivalent PRI particle. But significant differences in 
magnitude can also be expected for specific values of N r , simply because the Mie scattering 
coefficients a„ and b„ present distinct resonances for NRI and PRI particles [17]. 




Fig. 2. T x profile for an A-polarized TEMoo laser beam (X = 1064 nm, (U 0 = 3.7 um) and a spher- 
ical particle with a = 0.75 um and \N lm \ = \(T~ '. The position of the particle is (x 0 ,yo,Zo) = 
(0,1.9,0) um. Only the first 10 a„'s and b„'s significantly contribute to the torque profiles. Re- 
sonances in the Mie coefficients are reflected in the peaks observed. 
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Notice that, because of the ray optics characteristic alX «19.53, a significant number of 
Mie coefficients are necessary to account for a good description of the optical torque compo- 
nents. Although the MSCs can assume different complex values for PRI and NRI particles, it 
is only when a ~X or a « X that T x for a NRI sphere presents a distinct amplitude profile rela- 
tive to that of a PRI sphere. Figure 2 is a plot of T x as a function of N r for a = 750 nm, X = 
1064 nm, m 0 = 3.7 urn, N im = 10~ 7 , y 0 = 1.9 |im and z 0 = 0. Now, alX «0.705, so that the size 
parameter ka «4.43 and, consequently, only n up to 10 is used in (1) in order to ensure the 
adequate convergence of the GLMT [26,27]. Peaks at distinct \N r \ can be readily seen. 

Finally, Figs. 3(a) and 3(c) show three-dimensional views of T x as both N r and N im are 
changed, keeping the same beam parameters, y 0 and z 0 used in Fig. 2 (x 0 = 0). Except for some 
particular refractive indices, where resonances are observed, the optical torque first grows as 
Ni m increases from low values, thus reinforcing the fact that absorption is directly related to 
the SAM transfer. Then, over some specific ranges comprising intermediate Af,- m 's, T x resem- 
bles a plateau, decreasing after some limiting N im . This was originally observed for PRI par- 
ticles and, according to Fig. 3, such a behavior of the radial torques for linearly polarized 
TEMqo beams can also be extended for lossy NRI particles. In contrast, Figs. 3(b) and 3(c) 
reveal how the amplitude profile of T x is changed for specific liV r l where resonances may or 
may not occur. For example, t\ It\ is of the order of 10 4 over the range 10~ 8 < 

N im < 10~ 5 , which would mean a faster rotation if in both cases the particles had the same 
mass and geometry. Similar results for T y are observed if the particle is shifted along x. 

3.2. Circular Polarization 

In an optical tweezers system employing focused Gaussian beams for optically trapping bio- 
logical molecules and PRI particles in general, SAM is transferred from the incident photons 
to an absorbing particle located at the trap focus leading to a nonzero longitudinal torque T z 
[20]. This causes the trapped particle to rotate counter- or clockwise along z depending on the 
handedness of the polarization. 

Although the previous analysis for linear polarized light reveals different radial torque 
profiles and amplitudes for particular values of N r for PRI and NRI particles, it is hard to ob- 
serve such torques in a real experiment, basically because even NRI particles would neces- 
sarily be attracted towards the beam waist centre or repelled away from it [17]. In this way, T z 
would be the optical torque component responsible for making a trapped NRI or PRI particle 
rotates [20,23,30]. 

Let us again assume a first-order Davis description for a right-hand circularly polarized 
focused Gaussian beam with X = 384 nm, co 0 = 3.7 am propagating along + z. The particle is 
displaced along z and (xo,yo) = (0,0) with parameters a = 7.5 urn and N- m = 10~ 7 . Figures 4(a) 
and 4(b) show the optical torque T z and the radiation pressure cross-section along z, C pr>z , re- 
spectively, for five positive values of N r , the same adopted in Fig. 3 plus N r = 1.005. Radia- 
tion pressure forces along + z and -z are represented by the conditions C prz > 0 and C pr z < 0. 
The additional \N r \ = 1.005 was chosen because, in the GLMT, it is known to provide a region 
of negative C pr _ z (the undefined region -2.8xl0~ 4 m < z 0 < -1-OxlO" 4 m in Fig. 4(b)) and, 
therefore, a theoretical three-dimensional trap [37]. Because radial forces are null at any point 
along the optical axis, points where C WiZ = 0 in Fig. 3(b) represents the theoretical stable equi- 
librium points. For the parameters of Fig. 4, not a single particle with N r <0 would be trapped 
in a three-dimensional fashion, but yet a careful choice of the incident beam and the NRI par- 
ticle can eventually furnishes 3D-trappable NRI particles [17]. 
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Fig. 3. T x profile as function of N, m for (a) and (b): a PRI particle; (c) and (d) a NRI particle. 
The electromagnetic parameters of the beam and the location of particle are the same as in Fig. 
2. Different resonances in the Mie scattering coefficients lead to different torque profiles. 

From Figs. 4(a) and 2(c), one concludes that no significant differences in T z can be ob- 
served for the NRI and PRI curves using the adopted parameters, even though C pr . z may 
present profound changes in amplitude or in shape (Figs. 4(b) and 4(d)). Figure 5 is the equiv- 
alent of Fig. 4 for the particle of Fig. 2 located at (x Q ,y 0 ) = (0,0), where the effects of reson- 
ances in the Mie coefficients are pronounced and higher peak amplitudes in T z can occur. For 
N r = 1.005 (Fig. 5(b)), negative C pnz is found in the approximate range -6.5xl0~~ m < zo < 
-3.5xl0~ 5 m. Note from Figs. 4(d) and 5(d) that, for N r < 0, longitudinal radiation pressure 
cross section profiles do not vary significantly as the real part of the refractive index of the 
particle is changed, a fact already observed in a previous work [17]. Finally, we observe that, 
according to Figs. 3(c) and 3(d), the resonance on the MSCs increases the amplitude of the z- 
component of the torque for N im = 10~ 7 and N r = -1.480 or -1.071. This is why, in Fig. 5(c), 
the torque for these parameters is about two orders of magnitude more than that for N r = 
-1.19, for example, for which the geometry and electromagnetic properties of the scatterer 
does not reflect any resonance on the MSCs (for further details on these resonances on MSCs 
see Ref. [38]). 
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Fig. 4. r ; profile (due to circular polarization) for (a) PRI and (c) NRI particles with a = 7.5 
urn, Ni m = 1(T 7 and different N r . The associated radiation pressure cross sections C pr . : are shown 
in (b) and (d), respectively. In (b), the undefined region — 2.8xl(T 4 m < zo < -1.0x10 m for N, 
= 1 .005 represents negative C pl ■- not shown due to the logarithmic scale. Beam parameters are k 
= 384 nm and a> 0 = 3.7 um. 





Fig. 5. Same as Fig. 4, using the same parameters for both the laser beam and the particle as 
those of Fig. 2. 

Figure 6 is a three-dimensional view of T z as the complex refractive index of the particle is 
changed. The right-hand circularly polarized Gaussian beam has the same wavelength and 
beam waist as before, while the radius of the scatterer is fixed at a = 0.75 urn. The resonances 
at N r «- 1.480 and -1.071 for the NRI case are readily identified and, even at (x 0 ,y 0 ,Zo) = 
(0,0,0), a plot of T, versus \N, \ would result in a figure very similar to Fig. 2 for T x . 
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Fig. 6. 71 profile for an right-hand circularly polarized TEM 0 o laser beam (A = 1064 nm, w = 
3.7 um) and a spherical particle with a = 0.75 nm located at (;to,.Vo,Zo) = (0,0,0). Only the first 
10 a„ and b„ significantly contribute to the torque profiles. Resonances in the Mie coefficients 
are reflected in the peaks observed. 

The peak amplitudes in Fig. 2 and still valid for T- are essentially a consequence of a 
combination of the first ten a n 's and b„'s. Suppose, for example, aNRI particle with a = 0.75 
pm, N im = 10~ 7 and \N r \ «- 1.071. For this refractive index, the optical torque peak amplitude is 
basically due to a 7 , whereas the peak for \N r \ «-1.48 comes from b 9 , as shown in Fig. 7. In 
fact, we could say that these peaks arise from specific Mie scattering coefficients, all of them 
weighted by the beam shape coefficients of the incident beam. 




— 

— w 



-1.6 -1.5 -1.4 -1.3 -1.2 -1.1 




;<<=) 



— I", 
— - iw 



Fig. 7. (a) T : profile for an right-hand circularly polarized TEM 0 o laser beam (I = 1064 nm, w = 
3.7 um) and a NRI spherical particle with a = 0.75 nm located at (x 0 ,Vo,Zo) = (0,0,0). (b), (c) and 
(d) show the Mie scattering coefficients responsible for the peaks observed in (a). 



4. Conclusions 

The study of optical torques in optical trapping systems is extremely important and serves as 
an useful theoretical tool for predicting whether some biological particle will rotate, about 
some specific axis, under the presence of some arbitrary incident beam. In this way, the 
GLMT is an essential mathematical formulation to account for numerical optical torque calcu- 
lations because it can be used to describe the linear or angular momentum transfer from any 
laser beam to an arbitrary particle in any optical regime. The integral localized approximation 
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reduces computational time in the sense that it eliminates the undesirable and time consuming 
quadratures with double or triple integration. 

The Mie scattering coefficients present different phases and amplitudes depending on the 
geometry and the electromagnetic properties of the scatterer, and this is also true when we 
suppose an absorbing negative refractive index spherical particle. In this situation, new reson- 
ances appear which reflects our results for optical torques due to the polarization of the inci- 
dent beam. The inclusion of losses serves to make our model for the NRI particle more physi- 
cal, and is fundamental in SAM transfer analysis. But we do may expect that our previous 
analysis for lossless particles still applies for NRI spherical scatterers with very low losses 
(either due to the dispersive nature of these metamaterials or because of gain) at the operating 
frequency of the laser beam. 

The focused Gaussian beams explored here are not capable of transferring orbital angular 
momentum due to its azimuth symmetry, and we can naturally expect that other types of laser 
beams such as Laguerre-Gaussian and higher order Bessel beams, for example, will induce 
new optical torques in NRI particles. 

Experimental verification of our results is still a challenge because of actual technological 
limitations in nanofabricating effective homogeneous negative refractive index spherical par- 
ticles, especially in the optical regime. Although this may seem a little frustrating, it would be 
possible, in principle, to design a delicate experiment using macrostructures with a 2D NRI 
response in microwaves, small enough and with such a mass that, when impinged by a well- 
designed laser beam with sufficient power, it is mechanically oriented in a given plane as pre- 
dicted by our recent studies. 
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